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1. INTRODUCTION 

This document presents the data analysis procedures proposed for use with COMP- 
TEL instrument aboard the Compton Gamma Ray Observatory in the search for radia- 
tive neutrino decay from supernovae. The proposed analysis methodology is an extension 
of a standard procedure used by the COMPTEL team in searching for a variety of source 
types. We have applied the procedures to a set of simulated data to demonstrate the 
feasibility of the method to this project. 

2. v DECAY SIMULATION 


A computer model has been developed which simulates the radiative decay of a mas- 
sive neutrino species. Specifically, the simulation generates neutrino spectra relevant for 
typeTI supernovae and allows the neutrinos to stream from the surface of a proto-neutron 
star and decay in flight. An expanding spherically-symmetric shell is assumed and the 
appropriate relativistic kinematics are used. We have shown [1, 2] that the gamma ray 
spectra and intensities are dependent upon when in its history the supernova is observed, 
the supernova distance, the length of the observation period, and the postulated neutrino 
mass, lifetime and radiative decay branching ratio (m y ,r y ,5 7 ). 

Thus for a given COMPTEL supernova observation the first step in the data analy- 
sis chain is to generate gamma ray energy /angular distributions for a range of neutrino 
mass, lifetime and radiative decay branching ratios which correspond to a 3-D grid in the 
)-parameter space. These energy spectra and angular distributions can then 
be compared with the relevant COMPTEL data to search for the existence of radiative 
neutrino decay, and hence the existence of a massive neutrino species. 

3* MAPPING v SIMULATION TO COMPTEL DATASPACE 


Our simulation of radiative neutrino decay from supernovae produces a gamma ray 
energy spectrum and an angular distribution for given values of m Uy t v and B 1 . However, 
the comparison of the expected energy spectra from radiative neutrino decay with the 
COMPTEL data is non-trivial. As is well known the COMPTEL instrument does not 
measure the energy spectrum of a source directly. The relevant parameter for a Compton 
telescope is the scattering angle <j> which is defined through 
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where E\ and _E 2 is the energy deposition in the Dl and D2 modules, respectively. There- 
fore, to search for the decay gamma ray emission from radative decay we must map the 
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energy and and angular distributions into the COMPTEL dataspace. 

We take the approach described by de Boer, et al[3] . The COMPTEL data space 
consists of three parameters: x, ip, (j>. The first two parameters represent an aribtrary 
spatial coordinate system centered on the COMPTEL pointing axis and the third the 
scattering angle defined above. The COMPTEL dataspace response to a source (in the 
locally flat approximation) is a cone with apex at (Xo>^o)j the source position, running 
at an angle of 45° with the (x,t/>)-plane. However, because of finite resolution of event 
location within the detector modules, and energy measurement and deposition (which 
determine the scattering angle (j>) the cone becomes a mantle. The density of events 
along this mantle is given, for a particular gamma ray energy i? 7 , by the Klein-Nishina 
cross section. This mantle defines the PSF, /(x, <fi] Xo> V’o, E y ). 

The expected number of events due to a point source located at (x 0 , V*o) with intensity 
distribution (or energy spectrum) I(xo,i>o] £ 7 ) is given by 

e (x> V, V; Ei) = six, V) J i{xo , Vo; E y ) A(x a , Vo; £ 7 ) T fix, V, V; Xo, Vo, £ 7 ) dE 1 (2) 

while for a source with an angular distribution different from a point source 

e(x,V,V;£ 7 ) = ff(x,V) J Hx ,^',E y )A(x ,ip'\E^)T /(x,V,V;X ,V' ,E y )dx dip' dE y 

(3) 

where, g(x , VO is a geometrical absorption factor due to the finite size of the instrument 
(assumed to be independent of </>), A(x > 'V 1 Ey) is the effective area for a given photon 
arrival direction (x , )> an d T is the integration time. The discretized form of this is 

e(d) = g(d ) Y, Y fid, EJI( S] E y )X( S] E y ) (4) 

E 1 * 

where, d — (x^VsV) a dataspace bin, s ~ (x\V0 a sky pixel bin, and X = AT the 
exposure. Thus the second step in the data analysis chain is to map the output of the 
v decay simulation into the COMPTEL dataspace using the procedure outlined above. 
Once the expected number of events for each dataspace bin is computed the search for 
radiative neutrino decay can proceed as discussed below. 

4. SEARCHING FOR A SOURCE 

To compare our v decay simulation (after mapping into the COMPTEL dataspace) 
with the actual observational data we will use the Likelihood Ratio Method (LRM). Using 
our decay simulation and the mapping procedure described above we can determine a 
priori the number of gamma rays expected for each dataspace bin given a particular 
COMPTEL observation period, an assumed value of and B y (or range of values), 

and a source distance (redshift). In general the expected number of events per dataspace 
bin is given by (see section 3): 

,,, _ f B(d) if E 0 

1 ^ \ Bid) + g(d)^E 1 E a fid,s-,E^I(s > E y )X{s-E y ) if H v 
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where B(d) is the dataspace event density due to galactic diffuse, instrumental, and 
known source emissions, and H 0 and H v represent the hypotheses for background alone 
and background plus neutrino decay, respectively. With these assumptions we can assign 
a likelihood to a dataspace of n bins under a given hypothesis H as 

L({ni}] H) = JI e? exp (-ej/riil (5) 

1=1 

where n* is the observed number of counts in bin i . The likelihood ratio is defined by 

R = L{{n i } ] H„)/L{{n i } ] H 0 ) (6) 

and it has been established that A = 2 log R will adopt a x 2 distribution with 3 degrees 
of freedom (ie.,m^,7V,j9 7 ) in the event that the hypothesis H 0 is true. Thus the third 
step in the data analysis chain is to form the parameter A for each sampled grid point 
in the (m Ul T Uy B y )- -parameter space. 

5. ANALYSIS PROCEDURE 

For a given supernova, gamma ray energy spectra/angular distributions are generated 
for values of m„,7V,j0 7 based on the sourced distance (redshift) and a specific COMPTEL 
observation period. Using the methodology outlined above we can test for the presence 
of radiative neutrino decay from a specific supernova as follows: 

1. Determine the expected number of events per dataspace bin for hypothesis H 0 
(which includes background and any known sources). 

2. Compute the likelihood function L({n;}; H 0 ). 

3. Given a value of and 5 7 compute the expected number of events per datas- 

pace bin. This represents hypothesis H u , 

4. Compute the likelihood function i({n t }; H u ). 

5. Calculate A. 

6. Repeat steps 3 — » 5 for all values on the 3-D grid in R 7 )-parameter space. 

7. Find grid point with maximum A. This is the best fit point. 

8. If lambda max > A c (TBD) we reject the null hypothesis (H 0 ) y and can exclude 
grid points in (m U9 r Ui f? 7 )-parameter space are if they satisfy the condition A < 
A m ax “ X3' 

9. We claim a positive detection of radiative neutrino decay if A max > 5<r. 

By keeping a record of the mass and lifetime parameters sampled and the corresponding 
A values we can generate an exclusion plot in f? 7 )-parameter space for each super- 

nova observation, and as a final data product produce a exclusion plot for all supernovae 
observations combined. This method has the advantage that exposure and sensitivity 
differences between different observation periods, sources, and spacecraft status’ are in- 
cluded in the calculations and thus combining the results from multiple observations 
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Table 1: Sources/Observation Periods 


SN 

Obs. Period 

Source 
RA DEC 

GRO z-axis 
RA DEC 

Off-Axis Angle 
(deg.) 

1987A 

6 

083.86 -69.27 

091.28 -67.27 

3.0 

1987A 

17 

083.86 -69.27 

083.48 -72.27 

3.0 

1993J 

216 

147.86 69.30 

143.66 71.46 

4.7 


should be straight forward. 

6. SUGGESTED SOURCES 

The astrophysical sources which will be used to carry out this search are super- 
novae of type-II. We restrict our search to type-II supernovae since it has been shown 
experimentally [4, 5] that neutrinos are emitted during the collapse and cooling of the 
progenitor. In addition, the measured flux and energy of the electron- antineutrino burst 
from SN1987A was essentially the same as the predictions made with computer simu- 
lations of core-collapse supernovae. Thus we assume that, at least in neutrinos, type-II 
supernovae are standard candles. 

To sample the largest region of i? 7 )-parameter space the search should be 

performed with a combination of nearby and extra-galactic supernovae. Observations of 
SN1987A will allow a very sensitive search for massive neutrinos >few*100 eV because of 
its relatively close proximity, while extra-galactic supernovae which occur within COMP- 
TEL’s field of view will allow a search for a less massive neutrino. Observations of the 
recent type-II supernova SN1993J in M81, at approximately 3.2 Mpc, will allow a region 
of parameter space to be sampled which would not be senstively sampled by any other 
known source. We therefore recommend that the search for radiative neutrino decay 
concentrate on the sources and COMPTEL observation periods given in Table 1. It is 
also possible that a region of the parameter space not sampled by SN1987a or SN1993J 
will be accessible if an extragalactic supernova occurs serendipitously within the COMP- 
TEL field of view (near the pointing axis). 

7. APPLICATION 

As an example of the analysis procedures described above, we have processed two 
of the suggested sources to determine the general results which can be obtained. It 
should be stressed that the goal here is not to duplicate or estimate the results of a more 
detailed analysis to be completed in the near future, but rather to study the feasibility 
and overall features of the analysis technique. 

For this analysis a number of assumptions have been made: a) the effective area 
of COMPTEL is fixed at 30 cm 2 , b) energy and angular resolution is perfect, c) the 
geometrical absorption factor is fixed at 1, and d) the diffuse gamma-ray background is 
assumed uniform and no other sources are present in the COMPTEL FOV. During this 
visit we hope to reconcile these assumptions and, with the help and experience of the 
COMPTEL team, incorporate accurate models of the energy, off-axis angle, and position 
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Table 2: Input Parameters - v Decay 


SN 

Obs. Period 

Distance 

(Mpc) 

1987A 

1993J 

6 

216 

0.055 

3.2 


resolutions into the simulation and analysis procedures. 

The relevant input parameters to the radiative neutrino decay simulation (not given 
in Table 1) are given in Table 2 for two sources, SN1987A and SN1993J. For both 
analyses the neutrino output (i/ T -| -Tv, say) of the supernova in neutrinos was assumed 
to be 6.3 X 10 58 MeV with a Fermi-Dirac temperature of 8 MeV. Figure 1 shows the 
exclusion region in (m^ ^-parameter space from the SMM observations of SN1987A [6] 
(currently the best limits) assuming a branching ratio to a radiative decay mode of 1. 
Figures 2 and 3 show the exclusion regions obtained in the simulations for both SN1987A 
and SN1993J. 

8. QUESTIONS 

Question 1. We have developed software routines which perform the anal- 
ysis procedures described above - Does the COMPTEL team have similar 
software which can be used to compare our results? 

Question 2. Do the tables for the geometrical absorption factor (g), the 
effective area (A), and angular and position resolution exist for the range of 
energies and/or observational periods for which we are interested? 

Question 3. We have used a dataspace of the size 80°x80°x60° for 
Based on the experience of the COMPTEL team what are the suggested bin 
sizes in COMPTEL dataspace? 

Question 4. How can we get access to the background models or estimates 
made for each of the viewing periods of interest? 
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